Phase Diagram of a Model for Diluted Magnetic Semiconductors 
Beyond Mean-Field Approximations 
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A lattice spin-fermion model for diluted magnetic semiconductors (DMS) is investigated numer- 
ically, improving on previously used mean-field (MF) approximations. Curie temperatures are ob- 
tained varying the Mn-spin x and hole n densities, and the impurity-hole exchange J in units of the 
hopping amplitude t. Optimal values are found in the subtle intermediate regime between itinerant 
and localized carriers. Our main result is the behavior of the Curie temperature at large J/t, where 
a "clustered" state is observed and ferromagnetism is suppressed. Formal analogies between DMS 
and manganites are also discussed. 

PACS numbers: 75.50.Pp,75.10.Lp,75.30.Vn 
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Diluted magnetic semiconductors based on III-V com- 
pounds are attracting considerable attention due to their 
combination of magnetic and semiconducting proper- 
ties, that may lead to spintronic applications [Q, ^|. 
Gai-^Mn^As is the most studied of these compounds 
with a maximum Curie temperature Tc~110 K at low 
doping x, and with a carrier concentration p=(n/x)<l 
due to the presence of As antisite defects [0 or Mn inter- 
sticials H . It is widely believed that this ferromagnetism 
is carrier-induced, with holes introduced by doping me- 
diating the interaction between S=5/2 Mn-spins. This 
Zener mechanism operates in other materials as well [||. 

In spite of the excitement around 110 K DMS, room 
temperature ferromagnetism should be achieved for po- 
tential applications, with logic and memory operations in 
a single device. For this reason, a goal of the present ef- 
fort is to analyze the dependence of Tc on the parameters 
x, p, and J/t, helping in setting realistic expectations for 
DMS potential technological applications. This goal can 
only be achieved with good control over the many-body 
aspects of the problem, and for this purpose lattice Monte 
Carlo (MC) techniques are crucial, improving on previ- 
ously employed MF approximations. Our results lead to 
an optimistic view in this respect, since Tc is found to 
increase linearly with x up to x~0.25. 

Our effort builds upon previous important DMS theo- 
retical studies ||, |[ g, 0, g, |[ |l(J . However, to analyze 
whether Tc can be substantially increased from current 
values, techniques as generic as possible are necessary. 
In particular, both the strong interactions and disorder 
must be considered accurately, with computational stud- 
ies currently providing the best available tools. For these 
reasons, our work differs from previous approaches in im- 
portant qualitative aspects: (1) Some groups use a con- 
tinuum six-band description of DMS pf. (2) Other the- 
ories assume carriers strongly bounded to impurity sites 
and employ Hartree-Fock approximations. (3) Dy- 
namical MF theory (DMFT) may not capture the per- 
colative character of DMS, with a random impurity dis- 
tribution and cluster picture |8|, |) . (4) Other approaches 



use MF uniform states g], or introduce a reduced basis 
in simulations |lC| ] . While the previous work is important 
in describing current DMS materials, our goal is to es- 
tablish the phase diagram of a DMS model avoiding MF 
approximations. 

For the above mentioned reasons, here a generic MC 
study of a lattice spin-fermion model for DMS is reported. 
The Hamiltonian is 



H 



~t 4tCjV + J^Sl-ffi, 



(1) 
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where c ia creates a hole at site i with spin a, and the hole 
spin operator interacting antiferromagnetically with the 
localized Mn-spin Sj is 07 = c\ a a a .pcif3. The carrier can 
visit any site of the lattice (assumed cubic [|ll| ) . The in- 
teraction term is restricted to randomly selected sites, /, 
with a S—5/2 Mn-moment. These spins are here consid- 
ered classical with |S/|=1, as widely assumed, allowing 
for a MC simulation similarly as in other contexts [ pjj . 
Photocmission and band calculations locate J/t near 2 
0, [l(], [ll| Rather than fixing parameters to cur- 

rent DMS values, here J/t is unrestricted, while x and p 
mainly vary in the range [0,1]. Approximations include 
the neglect of on-site Coulomb repulsion, valid at small 
x and p where double occupancy is unlikely. In addition, 
Mn-oxides investigations Q show that an intermediate 
or large J/t plays a role analogous to a Hubbard U/t 
at any x 0]p5[. At low x, the probability of nearest- 
neighbors (NN) Mn-spins is also low (0.0625 at £=0.25), 
justifying the neglect of an antifcrromagnetic (AF) Mn- 
Mn coupling. The hole motion is described by a one- 
band tight-binding model, while a more realistic model 
should include many bands as well as spin-orbit interac- 
tion [^|. Despite this simplification, our study considers 
the underlying lattice, absolutely necessary for a qualita- 
tive understanding of the DMS phase diagram. 

The MC technique used here is as in Mn-oxides in- 
vestigations it includes the full Exact Diagonaliza- 
tion (ED) of the hole sector for each MC spin con- 
figuration, and density-of-states expansion calculations 
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FIG. 1: (a) S(q=Q) vs. temperature, T, for 4 3 (circles), 6 3 
(squares), and 8 3 (diamonds) clusters at J/t— 2.0, x=0.1 and 
p=0.4, using the MC technique. The inset shows the magneti- 
zation \M\ for the same clusters, with a vertical scale referred 
to the asymptotic 8 3 T=oo value, (b) \M\ (circles) and spin- 
spin correlation C at maximum distance (squares) vs. T, on a 
4 3 cluster, at J/t=1.0, £=0.25, and p=0A. The 0.25 horizon- 
tal line indicates the \M\ T=oo asymptotic value, (c) S(q—0) 
vs. T, at J/t=2.0, £=0.25, and p=0.4 using 4 3 (circles) and 
6 3 (squares) clusters. In all cases Tc is indicated. 

just . The latter allows us to reach clusters with up 
to 8 3 =512 sites if up to 40 terms are included, reach- 
ing an accuracy comparable to ED for smaller clusters. 
Both methods are nearly exact, and the error bars of 
our results mainly arise from intrinsic thermal fluctua- 
tions and averages over several random Mn-disorder con- 
figurations. Comparing estimations of different clusters 
and based on previous experience with similar models 
PL Tc can be calculated within a ~25% accuracy, suf- 
ficient for our purposes [ fl7|| . The order parameter for 
the ferromagnetic-paramagnetic transition was taken to 
be the absolute value of the magnetization of the Mn- 

spins normalized to 1, namely |Af|=^ r(Si ■ Sr). 
Size effects are better visualized in the zero-momentum 
spin structure factor S(q=0)=-^ J^i r(Si-Sr)- Another 
useful quantity is the spin-spin correlation at distance d, 
C ( d )=W(d)l>2\i-R\=d(Si ■ Sr}> wh ere N(d) is the num- 
ber of pairs of Mn moments separated by a distance d. 

Typical results for small and intermediate J/t of our 
large-scale computational effort are in Fig. la. There 
S(q=0) and \M\ vs. temperature T are shown for three 
cluster sizes, x=0.1, and p—OA. Note the small size de- 
pendence of the magnetization (inset), and the volume 
growth of S(q=0) at fixed T<T C . The estimated T c /t 
is ~0.04, with an uncertainty 0.01 sufficiently small for 
our purposes. Even with just the 4 3 cluster, Tc could 
be estimated fairly well, as shown in Fig. lb. This is im- 
portant to simplify our computational search for optimal 
Tc's varying many parameters. In Fig. lb, the tempera- 
ture where a deviation from the high-T limit is found is 
slightly larger than the Tc/i=0.04 obtained from larger 



clusters (indicated). Studying the spin-spin correlation 
at the largest available distance, a nonzero value charac- 
teristic of an ordered ferromagnetic (FM) state was ob- 
tained at T just below 0.04. Figure lc provides another 
example of our comprehensive Tc study, using just two 
cluster sizes at the x-p location of our most optimal Tc, 
at fixed J/t=2. Here the use of only 4 3 and 6 3 clusters 
provides once again a fairly accurate value Tc^0.08t. 

To understand the qualitative Tc trends, first consider 
the simplest case: the p dependence at fixed J and x. Us- 
ing the results in Fig. 2a contrasted against Fig. lb (same 
cluster size) Tc is found to change by a factor ~2, when 
p varies from 0.1 to 0.4. However, this tendency does 
not continue with increasing p, since at p=l or beyond, a 
FM state is not formed: the Pauli principle reduces dras- 
tically the carrier kinetic energy, leading instead to an AF 
state. An example at p=3 and on an 8 x 8 cluster (results 
are qualitatively similar in two and three dimensions) is 
in Fig. 2b, where the oscillations in the spin correlations 
indicate staggered order. In general, the optimal p is 
~0.5, between the hole empty p—0 and saturated p=l 
limits, as found with DMFT [Q. A similar result occurs 
in Mn-oxide models, recovered from Eq.(l) at x=l. In 
that context, investigations at large Hund coupling, the 
analog of J for DMS, have shown that p=0.5 optimizes 
Tc to a number ^0.11-0.13t |8), likely an upper-bound 
on the Tc that could be achieved with Eq.(l). 
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FIG. 2: (a) Magnetization \M\ (circles), and spin-spin corre- 
lation at maximum distance C (squares) vs. T, at J/i=1.0, 
2=0.25, and p=0.1, using a 6 3 cluster, (b) C(d) vs. d at 
p~3.0, 2=0.25, J/t=1.0, and T/i=0.005, using an 8 2 cluster. 
The oscillations indicate an AF state, (c) \M\ (circles), spin- 
spin correlation at minimum distance C(dmin) (squares), and 
at maximum distance C(d max ) (diamonds) vs. J/t, for a 4 3 
cluster at 2=0.25, p=0.1, and T/t=0.005. (d) A£= £(0=tt)- 
E(0=0) vs. J/t calculated exactly on finite but large clusters 
at T=0 for a system of 2 Mn-spins and 1 electron, 9 being the 
relative angle between the Mn-spins. Results in 1, 2, and 3 
dimensions are indicated. The spin distance is such that the 
associated effective 2~0.1 is the same in all cases. 
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Consider now the J/t dependence of Tc- The MF ap- 
proximation suggests T™^ocJ 2 . However, this does not 
hold when more accurate methods are used in the cal- 
culations. In fact, for J/t— >oo and a Mn dilute system, 
the holes are trapped in Mn-sites, reducing drastically 
the conductance and Tq- Small FM clusters of spins 
are formed at a temperature scale T*, but there is no 
correlation between them, leading to a global vanish- 
ing magnetization ||. These results cannot be obtained 
within a mean-field approximation. The large- J/t ideas 
can be tested in our MC simulation by monitoring the 
short- and long-distance behavior of the spin-spin corre- 
lations C(d). In a "clustered" state (large J/t), C{d) at 
the shortest distance can be robust at T<T*, but C(d) 
at the largest distance vanishes due to the uncorrelated 
nature of the magnetism between independent clusters 
(see Fig. 2c). This subtle effect explains the incorrect MF 
prediction, since T™'' ~T*, which grows with J/t, rather 
than the true Tc (see also Fig. 4a). Since both in the 
J/i~0 and J/t=oo limits Tc is suppressed, an optimal 
J/t\ op t must exist where Tq is maximized. Simulation 
results as in Fig. 2c indicate that the optimal J/t value 
is close to 2. This phenomenon is not captured in itiner- 
ant g] or localized || limits nor by DMFT @, but it is 
observed in the present generic MC simulations. 

The existence of a J/t\ opt can be illustrated just using 
two spins and one carrier in a finite cluster at T—0. For 
any fixed angle 9 between the Mn-spins, assumed copla- 
nar, the energy is found exactly. The ground state of 
this p=0.5 system is always at 9=0 (FM-), while the en- 
ergetically worse state is 9=ir (AF-configuration) . Their 
energy difference AE is a crude estimation of the FM 
state stability (Fig. 2d). An optimal J/t is found in all 
dimensions, with stability increasing with the coordina- 
tion number PQ| . The result Fig. 2d is understood mea- 
suring the electronic density n(i) of the same problem on 
a chain (Fig. 3a). At small J/t, the derealization mani- 
fests in the nearly uniform density, leading to weak FM. 
At large J/t, strong localization decouples the Mn-spins, 
producing again weak FM. However, there is an optimal 
value where the system takes advantage of J/t, but also 
allows for a nonzero effective coupling among separated 
classical spins, leading to a stronger FM. 

Consider now the x dependence of Tc- For simplicity, 
J/t=2 is mainly studied, which is both close to optimal 
and experimentally realistic p2| . Fig. 3b shows Tc vs. x 
at p~0.4, and for two reasonable values oft. Experiments 
|jj indicate a linear growth of Tc up to 5% (shown), as 
in the numerical results. The slope of Tc vs. x is in 
remarkable agreement with MC predictions, in a reason- 
able range of t. Regarding x>0.05, a reduction of Tc was 
originally reported in experiments 0]. However, recent 
data gathered with an optimized annealing treatment f23f| 
indicate a Tc "plateau" . This still seems in contradic- 
tion with the linearly growing Tc of the MC results, but 
it suggests that even more refined thin-films may con- 
tinue increasing Tc with increasing x. The MC results 
clearly indicate linear behavior up to x~0.25 (Fig. 3b). To 
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FIG. 3: (a) Exact T=0 local carrier density, for 1 carrier 
and 2 parallel spins at sites 10 and 20, of a 30-site chain, 
varying J/t. (b) Tc vs. x for Gai-^Mn^ As [Ell (thick line), 
and for model Eq.(l) using t=0.3 eV (squares) and t=0.5 eV 
(circles). In both cases, p~0.4 and J/t=2.0. Typical error 
bars are shown. Stars are results for p^0.5, J/t=l, and t=0.5 
eV. (c) Magnetization \M\ vs. x for model Eq.(l) (circles) at 
T/t=0.005, J/t=2, p=0.4, compared with the experimental 
value (lj for Gai-xMnxAs at x~ 0.07 and T=2K (squares). 



the extend that our model describes DMS quantitatively, 
higher values of Tc could be expected experimentally. 
Regarding the presence of a Tc maximum at x=0.25: the 
origin of this effect is the growing probability with x of 
having both holes and Mn-spins at NN-sites. In this case, 
AF links are formed since J / 1=2 is not so strong to keep 
the link FM, reducing Tc at large x even at p=0.5. As 
J/t grows, the effect diminishes and the maximum in Tc 
moves toward x=l, as naively expected. Reciprocally, as 
J/t decreases from 2 (J/t=l shown in Fig. 3b), the maxi- 
mum in Tc moves toward smaller x's, and only t>0.5 can 
provide high-T ferromagnetism. This illustrates the key 
role that the optimization of J/t plays in these models, 
effect not captured by MF approximations. 

\M\ at T~0 is in Fig. 3c. In agreement with experi- 
ments, the x^O.l result indicates a magnetization ~50% 
of its maximum value. This nonsaturated behavior orig- 
inates in the random distribution of Mn-spins, since Mn- 
clusters are formed providing a trap to holes. Non- 
clustered spins are not much visited by those holes, and 
their spins are not polarized. With growing x, holes are 
more itinerant, polarizing the entire sample (Fig. 3c) p4| . 

In summary, MC investigations of a spin-fermion 
model for DMS unveils substantial differences with pre- 
viously reported results employing MF techniques. The 
subtle regime of intermediate J/t appears the most rel- 
evant in these compounds. Tc~0.08i is an upper limit 
for the FM critical temperature, result close to those ac- 
cepted for x=l §, ||. Our main results are summarized 
in Fig. 4, that contain (a,b) the nontrivial J/t dependence 
of Tc and T* , and (c) Tc with varying x and p, at op- 
timal J/t. To the extend that the present model is ap- 
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FIG. 4: (a) MC phase diagram in 2D varying J/t, at fixed 
x and p. At large J/t, a broad scale T* corresponds to the 
formation of uncorrelated clusters. Tc is the "true" transition 
temperature, defined as the T where FM correlations develop 
at the largest distance in the clusters used. At small J/t, those 
temperatures are similar. The optimal J/t is intermediate 
between itinerant and localized regimes, (b) Schematic phase 
diagram believed to be valid both in 2D and 3D, with the 
clustered and FM states indicated, (c) Numerically obtained 
Tc vs. x and p, at J/t— 2.0. Filled circles are from model 
Eq.(l) with t=0.5 eV. The open square corresponds to the 
experimental value for Gai-^Mn^As at x~0.1. 



plicable to DMS materials, broad guidelines to improve 
T c can be established: (i) The optimal J/t| opt ~2 must 
be intermediate between the itinerant and localized lim- 
its (Fig. 4a, b). This J/t, or larger, is expected to keep 
the semiconducting nature of the state at T>Tc- Only 
band calculations beyond our model can predict which 
particular material will have such an optimal J/t. (ii) 
x should be increased beyond 0.1. At J/t\ op t, the best 
value is x~0.25. Currently, x=0.14 is the experimen- 
tal limit (^6). (iii) The number of antisite defects must 
be controlled such that p^0.5 (p~l would be detrimen- 
tal due to competing antiferromagnetism). (iv) As the 
coordination number grows, Tc grows, (v) Finally, the 
simplest procedure to increase Tc relies on increasing the 
scale t. In fact, (Ga,Mn)As and (In,Mn)As have differ- 
ent hybridization strengths p7| , and this should be an 
important consideration in studying new materials. Our 
work also suggests formal analogies between DMS and 
manganite models, with similar Tc's, and a related clus- 
tered state above ordering temperatures. 
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Note also that the nonsaturated low-T state can be fully 
polarized with relatively small magnetic fields: in our 
studies, fields of just 10 T are sufficient to raise \M\ from 
0.5 to nearly 1.0, at a;=0.1. 

Lao.6Sro.4Mn03 has a Tc~340 K, and our simulations 
suggest that similar values could b e reached for DMS . 
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